Pulling absorbing and collapsing polymers from a surface 
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A self-interacting polymer with one end attached to a sticky surface has been studied by means of a 
flat-histogram stochastic growth algorithm known as FlatPERM. We examined the four-dimensional 
parameter space of the number of monomers up to 91, self-attraction, surface attraction and force 
applied to an end of the polymer. Using this powerful algorithm the complete parameter space of 
interactions and force has been considered. Recently it has been conjectured that a hierarchy of 
states appears at low temperature/poor solvent conditions where a polymer exists in a finite number 
of layers close to a surface. We find re-entrant behaviour from a stretched phase into these layering 
phases when an appropriate force is applied to the polymer. We also find that, contrary to what may 
be expected, the polymer desorbs from the surface when a sufficiently strong critical force is applied 
and does not transcend through either a series of de-layering transitions or monomer-by-monomer 
transitions. 

PACS numbers: 05.50.+q, 05.70.fh, 61.41. +e 



New experimental methods in the physics of macro- 
molecules [l| have been used to study and manipulate 
single molecules and their interactions. These methods 
make a contribution to our understanding of such phe- 
nomena as protein folding or DNA un-zipping; one can 
push or pull a single molecule and watch how it responds. 
It is possible to apply (and measure) forces large enough 
to induce structural deformation of single molecules. One 
can monitor the mechanism of some force-driven phase 
transition occurring at the level of a single molecule. The- 
oretical understanding of this behaviour has attracted 
much attention 0, 0, l1 

The response of a single polymer to an external force 
under good solvent conditions was considered some time 
ago 5]. The response under poor solvent conditions (be- 
low the #-point), where the self-attraction and an exter- 
nal force compete with each other, was examined later 
[lH H HH Another phenomenon commonly studied in 
polymer physics is the adsorption of a polymer tethered 
to a "sticky" wall. The response of such a polymer to a 
force perpendicular to the wall has also recently been con- 
sidered 0,EJE1- However, when both the self-attraction 
(ie. monomer-monomer attraction that leads to polymer 
collapse) and the surface attraction (ie monomer-wall at- 
traction that leads to adsorption) compete the response 
to an external force has not yet been elucidated (some 
interesting results can be found in [ll|). Certainly, the 
full phase diagram has not been considered. Making such 
a study now is all the more timely because of the very re- 
cent discovery [f3 | of a new low temperature phenomenon 
of layering transitions (without a force) . It is this layer- 
ing phenomenon that raises the intriguing question about 
the response a low-temperature polymer may have to an 
external force. In the layering state a polymer is tightly 



confined within a fixed number of layers above the wall. 
It may be therefore be especially interesting experimen- 
tally to examine such a situation. 

We demonstrate for the first time how the full two- 
dimensional phase diagram of surface and self-attraction 
changes as the force is increased. The desorbed extended 
regime, which changes its scaling behaviour as soon as 
the force is made non-zero, simply grows as the force 
is increased. The second-order phase transitions of ad- 
sorption and collapse become first order. The rest of the 
phase diagram remains relatively unaffected as long as 
the force is small. After the force passes a critical value, 
which depends on the zero temperature force required to 
pull a polymer from a wall, a re-entrant behaviour occurs 
at low temperatures. For different values of the force, this 
re-entrant behaviour occurs for both the adsorption and 
collapse of polymers, including the layering phases men- 
tioned above. We provide a full force-temperature dia- 
gram for all ratios of surface attraction to self-attraction. 
All this is achieved with the use of a recently developed 
algorithm, FlatPERM ^jjj], that is specifically designed 
to obtain information about the whole phase diagram in 
one simulation run: it is effectively a stochastic enumer- 
ation algorithm that estimates the complete density of 
states. 

The model we have considered is a self-avoiding walk 
on a three-dimensional cubic lattice in a half-space. The 
self-avoiding walk is attached at one end to the wall 
with surface energy per monomer of e s for visits to the 
wall. The self-avoiding walk self-interacts via a nearest- 
neighbour energy of attraction per monomer-monomer 
contact. (Note that the attractive energies Eb and e s are 
taken to be positive.) A force / is applied in the direction 
perpendicular to the boundary of the half-space (wall). 
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The total energy of a configuration tp n of length (number 
of monomers) n is given by 

E n {ip n ) = -m b (ip n )e b - m s (ip n )e s - fh(tp n ) (1) 

and depends on the number of non-consecutive nearest- 
neighbour pairs (contacts) along the walk mb, the num- 
ber of visits to the planar surface m s , and the height 
h in the direction perpendicular to the boundary (wall) 
of the half-space. Figure 1 shows a diagram of the 
two-dimensional analogue. For convenience, we define 
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FIG. 1: A diagram showing the two-dimensional version of 
the three-dimensional model simulated. 
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(3 b = e b /k B T, fa = e s /k B T and = f/k B T where T is 
the temperature and k B is the Boltzmann constant. The 
partition function is given by 

Z n ([3 b ,p s ,(3 f )= Cn, mb , me , h e 0bmb+p ° m ° +p f h (2) 

with C n ^ mhl m s ,h being the density of states. It is this den- 
sity of states that is estimated directly by the FlatPERM 
simulation. Our algorithm grows a walk monomer-by- 
monomer starting on the surface. We obtained data for 
each value of n up to n max — 91, and all permissible 
values of mb, m s , and h. 

When / = the phase diagram of the model contains 
several phases and transitions between them 0, 0, 0, 
Hif . For small /3& and (3 S there is a desorbed extended 
(DE) phase with the polymer behaving as a free flexible 
polymer in solution (ie. swollen or extended in three di- 
mensions). For /3b fixed and small, increasing (3 S leads to 
a second-order phase transition (adsorption) to a state in 
which the polymer is adsorbed onto the wall and behaves 
in a swollen two-dimensional fashion (AE). Alternately, 
if (3b is increased at small (3 S a second-order collapse tran- 
sition occurs to a state resembling a dense liquid drop. 
This phase is known as desorbed collapsed (DC) on the 
assumption that it has little contact with the wall [lilllfij . 
However, it has been subsequently argued that for 
larger (3b there is also a polymer-surface transition to a 




FIG. 2: Plots of the logarithm of the maximum eigenvalue 
of matrix of second derivatives of the free energy for three 
values of (3f = 0.0, 1.5 and 3.0. In the top far corner of the 
plot is the location of the desorbed- extended phase when (3f 
and the 'Stretched' phase when (3j > 0. The location of 
the 2-layer (2L), 3-layer (3L), adsorbed-collapsed (AC) and 
ads orbed- extended AE phases do not seems to move greatly as 
f3f is increased. 



Surface-Attached Globule (SAG) phase, where the poly- 
mer behave as a liquid drop partially wetting the wall. 
This transition will not be seen directly by studying ther- 
modynamic polymer quantities as it occurs as a singular- 
ity in the surface free energy and not the bulk free energy 
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of the polymer. When j3 a is large, so that the polymer 
is adsorbed onto the wall, increasing /3& will result in a 
two-dimensional (second-order) transition to a adsorbed 
and collapsed phase (AC). In recent work 0] this AC 
phase was also referred to as the 1-layer phase because 
for very large (3b and (3 S < (3b there exist meta-stablc 
Mayer phases where the polymer is two-dimensionally 
collapsed and more-or-less restricted to t layers parallel 
to the wall (for small £) . A series of first-order transitions 
between adjacent values of £ occur as (3 S is varied at fixed 
(3b- All these transition lines can be seen in the Figure 2 
(a) which shows a plot of the logarithm of the maximum 
eigenvalue of the (2 x 2) matrix of second derivatives in 
the variables f3b and (3 S of \og{Z n ([3b, (3 S , /?/)) for fixed 
f3f = 0. The local maxima indicate transitions. The 
transition to the AC phase is expected to be first order 
in the thermodynamic limit. 

Using the evidence available in the literature 0, EJ 0, 
EJ E3 , let us now consider what we can expect when 
/ > 0. The first important feature to note is that the 
isotropic DE phase is replaced by an anisotropic phase in 
which the height of the end point of the polymer scales 
linearly with n; we denote this phase as the Stretched 
phase. Consequently the transition from stretched to ad- 
sorbed phases becomes first-order 0. Likewise, at least 
in three dimensions Q , the transition from the vertically 
stretched phase to the collapsed phase also becomes first- 
order. This implies that the multicritical point (where for 
/ = the DE, AE and DC phases meet) is now a triple 
point: the meeting of three first order lines. The transi- 
tion from the AE to AC phases should not be effected by 
the application of a small force as it acts in a direction 
perpendicular to the driving phenomenon of planar col- 
lapse. Finally, it is intriguing to ask what happens to the 




FIG. 3: A typical configuration resulting from the application 
of the critical force f c to a polymer in the 2-layer adsorbed 
collapsed phase. 



layering phases observed in [12j. One can imagine that 



the force simply extends a vertical 'tail' from a layered 
block (see figure 3) and that as the force is increased the 
monomers are peeled off one at a time with correspond- 
ing micro-transitions fll| for each monomer pulled until 
a vertical rod is achieved. Instead we see at some point a 
sharp first order transition between the highly stretched 
vertical rod and a layered system with short tail. 

In figures 2(b) and 2(c) we show plots of the 
logarithm of the maximum eigenvalue of the matrix 
of second derivatives in the variables (3b and (3 S of 
\og(Z n ((3b, (3 S , (3f)) at fixed (3f (as in Figure 2(a)) but at 
values of (3f being 1.5 and 3.0. It is clear that as (3f is in- 
creased the stretched phase that occurs for small (3b and 
(3 S expands while the positions of the other phases and 
transitions move little. We immediately note that these 
plots do not tell the whole story since physically one is 
usually interested in fixing the force / rather than (3f. fix- 
ing (3f implies that the force applied goes to zero at low 
temperatures. It is for this reason that the re-entrant be- 
haviour for absorbing polymers 0, Q is not seen directly 
in these plots. However, re-entrant behaviour does occur 
and occurs for any ratio of surface to bulk energies. Let 
us now consider the more traditional force-temperature 
diagram and return to this point. 

In figure 4 a plot of the force f c (T, a) needed to pull a 
polymer from the wall for a ranges of temperatures and a 
parameter a which measures the relative strength of the 
surface (wall) interaction to the self-interaction. We have 




FIG. 4: A plot of the force f c needed to pull a polymer from 
the surface against temperature T and a parameter a. The 
parameter a controls the relative strength of wall attraction 
and self-attraction with e s ~ a and e\, = 1 — a. The limiting 
cases of surface desorption and of pulling a collapsed polymer 
are easily visible in the plot for a = 1 and a — 0, respectively. 

parameterised the energies of surface and self-attraction 
as e s = a and eb = 1 — ot respectively. Using this pa- 
rameterisation for < a < 1 gives the whole range of 
attractive activities: the ratio of surface to bulk activ- 
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ities is given as /3 s //3& = a/(l — a) and so is constant 
for fixed a. For a = we have e s = and Eb = 1 
which corresponds to pure self-attraction while the other 
boundary of the parameter space with a = 1 gives e s = 1 
and £b = which is the pure surface adsorption case. 
This extends the diagrams given in 0, in which only 
adsorption is considered. If a force smaller than f c is 
applied the polymer is in the phase appropriate to the 
value of a: either collapsed or adsorbed or both. On 
the other hand for forces larger than f c the polymer 
is in the 'stretched' phase. We immediately note that 
the reentrant behaviour observed in the adsorption-only 
case 0, persists for all a. Fixing the force to be at 
a value slightly larger than the zero temperature criti- 
cal force and then increasing the temperature leads to 
transitions from the stretched state to a non-stretched 
phase and back again to the stretched state. This arises 
due to the entropy of the zero temperature state; one 
can easily extend the arguments in to demonstrate 
that re-entrant behaviour can occur even when the zero- 
temperature configuration of the non-stretched state is a 
Hamiltonian (fully compact) walk in a cube rather than 
a totally adsorbed polymer. 

Let us discuss the re-entrant behaviour returning to 
the (/3b,/3 s ,/3/)-parametrisation. For fixed Eb, e s , and /, 
changing the temperature T implies moving on a ray in 
the (/3&,/3 s ,/3/)-space. At high temperatures the system 
is in the stretched phase near the origin. At low tem- 
peratures, the state of the system depends on the choice 
of £t,, e s , and /. For very large / the system remains 
stretched at all temperatures, whereas for very small / 
and low temperatures, the system is in a layered phase. 
The re-entrant behaviour manifests itself in the following 
way: there is a range of / for which the system, when 
moving along a ray in the /3 s ,/3/)-space (ie. for fixed 
Eb, e s , and /), changes from a stretched state near the 
origin to a layered one at intermediate temperatures and 
then changes back to a stretched state as 0b, f3 s , and /3/ 
become larger. 

If the critical force is zero then the curve in the T — a 
plane corresponds to the phase boundary of the DE phase 
with the apex of the curve around a ps 1/2 being the lo- 
cation of the multicritical point. On the other hand for 
T = there is a kink in the function f c (a) at exactly 
a = 1/2 which is a consequence of the first order point 
coming from the transition from SAG/layer phases from 
small a to the AC phase for larger a. There is the ap- 
pearance of a kink joining the multicritical point to the 
zero temperature transition which is presumably a finite 
temperature effect of the transition to the AC phase. 

In this paper we have studied how the phase diagram 
of a self-attracting polymer that is also attracted and 
tethered to a flat wall changes as a vertical force is ap- 
plied to the un-tethered end of the polymer. We have 
accomplished this using a flat histogram Monte Carlo 
simulation that is capable of studying the whole range 



of microscopic energies, temperature and polymer length 
up to a maximum of 91 monomers. We demonstrate that 
re-entrant behaviour occurs at low temperature and for a 
range of forces for all relative strengths of self-attraction 
and surface attraction. For small forces we have found 
that only the transition boundary of the "stretched" 
phase moves with increasing force while the rest of the 
phase diagram is relatively unchanged. In contradiction 
to what may be expected we have found that the novel 
layering meta-phases found for large but finite polymer 
length are unaffected by small forces. 
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